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ABSTRACT 

We use high resolution N-body/SPH simulations to study the hydrodynamical inter- 
action between the Large Magellanic Cloud (LMC) and the hot halo of the Milky 
Way. We investigate whether ram-pressure acting on the satellite's ISM can explain 
the peculiarities observed in the HI distribution and the location of the recent star 
formation activity. 

Due to the present nearly edge-on orientation of the disk with respect to the orbital 
motion, compression at the leading edge can explain the high density region observed 
in HI at the south-east border. In the case of a face-on disk (according to Mastropietro 
et al. 2008 the LMC was moving almost face-on before the last perigalactic passage), 
ram-pressure directed perpendicularly to the disk produces a clumpy structure charac- 
terized by voids and high density filaments that resemble those observed by the Parkes 
HI survey. As a consequence of the very recent edge-on motion, the Ha emission is 
mainly concentrated on the eastern side where 30 Doradus and most of the supergiant 
shells are located, although some Ha complexes form a patchy distribution on the 
entire disk. In this scenario only the youngest stellar complexes show a progression in 
age along the leading border of the disk. 
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1 INTRODUCTION 

The Large Magellanic Cloud (LMC) has revealed a very 
complex structure both in the stellar and in the gaseous 
component. The elongation of the stellar disk in the direc- 
tion of the Galactic center, its substantial vertical thickness, 
the warp and the strong asymmetric bar are naturally pre- 
dicted by numerical simulations as a result of the gravita- 
tional interaction between the LMC and the Galaxy (Bekki 
& Chiba 2005, Mastropietro et al. 2005, hereafter M05). The 
old stellar distribution appears to be quite smooth in the 
outer parts of the disk, with no signs of spir al structures 
out to a radius of 10 kpc (van der Marell lioOll V Within the 
same radius the HI large scale structure reveals the pres- 
ence of several asymmetric features that have no equivalent 
in the old stellar disk. The gaseous disk is characterized by 
the presence of an elongated region located at the south- 
east of the galaxy and aligned with the border of the optical 
disk, whe re the column density distr i bution shows a stee p 
increase (jStavelev-Smith et al.1 l2003l : IPutman et aDl2003l l. 



* E-mail: chiara.mastropietro@obspm.fr 



Since the LMC proper motion vector is directed to the east, 
it appears natural to associate this high density region with 
ram-pressure acting on the leading edge of the disk due to 
the orbital motion of the LMC and its consequent interac- 
tion with the diffuse hot gas in the halo of the Milky Way 
(MW). 

The presence of an extended hot halo surrounding 
galaxies and in hydrostatic equilibrium within the dark 
matter potential is expected by current models of hierar- 
chical structure formation (Mail er Bull ock 2004). In the 
MW, X-ray absorption lines produced by hot (T ^ 10^ 
K) gas are detected i n the spectra of s everal bright AGN 
(jWilliams et al.1 l2005l : iFang et al.l l2006h . Some ionization 
features discovered in the Magellanic Stream and high veloc- 
ity Clouds indicate that this distribution of hot gas extends 
well beyond the Galactic disk ( > 70 kpc). Constraints from 
dynamical and thermal arguments fix its density in a range 
between 10~^ and 10~^ cm ~^ at the LMC distance from the 
Galactic center (but Kaufmann et al. 2009 suggest a value 
ten times higher). 

ICioni et al.l (|2006[ ) have performed a detailed analysis of 
the LMC global star formation rate using asymptotic giant 
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branch stars. They find an irregular and patchy distribution 
in age, with the youngest carbon-rich systems located at the 
south-east of the disk. The present star formation activity is 
rather clumpy and concentrated in stellar complexes char- 
acterized by intense HII emission and associated with bright 
Ha filamentary bubbles. Most of these very young structures 
lie on the south-east of the disk, in the proximity of 30 Do- 
radus, the largest star forming region of the LMC, some are 
located in the bar and the remainder form an asymmetric 
pattern that covers the entire disk with no apparent relation 
to the global geometry of the satellite. 

It is not clear which is the overall physical mechanism 
responsible for triggering star formation with the observed 
asymmetric pattern and different models have been pro- 
posed in the past. The s tochastic self-propagati ng star for- 
mation (SSPSF) model (jCerola k Seidenlll978h predicts a 
clear age gradient in the LMC's stellar complexes, with the 
edges being younger with respect to the ce n ter, in contra- 
dictio n with observations (|Braun et al.lll997l l . Ide Boer et al.l 
proposed a scenario where the bow shock originated 
by the motion of the LMC through the hot galactic halo 
compresses the leading edge of the disk and induces star for- 
mation. The pressure at the south-eastern edge of the LMC 
is ind eed 10 times higher th an the average in the rest of the 
LMC (jBlondiau et al.lll997l ). This model, which assumes the 
orbital motion vector lying in the plane of the disk, predicts 
increasing ages of the stellar complexes in the direction of 
the rotation, due to the fact that the material compressed at 
the front side of the disk moves, in time, away to the side. 
The youngest systems would indeed lie at the south-east 
border of the disk, where the relative velocity between the 
corotating interstellar medium and the external diffuse gas 
is maximum. Several giant structures along the outer east 
and north edge of the LMC actually show a progression in 
age in a clockwise direction: moving from south-east to the 
north LMC 2, 30 Doradus and LMC 3, LMC 4, NGC1818. 
In particular the difference in age between 30 Doradus and 
LMC4 is exactly their distance along the border of the disk 
divided by the satellite's rotational velocity (Harris, private 
communication). Grebel & Brandner (1998) studied the re- 
cent star formation history of the LMC using Cepheids and 
other supergiant stars and found that although the majority 
of the star formation events in the last 30 Myr are concen- 
trated on the east border, others are distributed across the 
entire disk in partial contrast with the bow shock induced 
star formation model, that can not explain them. 

In this work we use high resolution SPH simulations to 
study the effects of the interaction between the LMC inter- 
stellar medium and the diffuse hot halo of the MW. We want 
to investigate whether the ram-pressure acting on the lead- 
ing edge of the LMC disk is responsible for the increase in 
density observed in the south-east and for triggering star for- 
mation. The analytic model of de Boer et al. (1998) assumes 
a pu re edge-on model, but according to Ivan der Mare l et al.l 
(|2002. ) the present angle between the LMC disk and the 
orbital motion is nearly 30°. Even in the absence of preces- 
sion and nutation, this angle is subjected to large variations 
during the orbital period in such a way that compression 
produced by the external hot gas can affect in time both 
edge-on and face-on. 

Moreover, the ram-pressure felt by the LMC is not con- 
stant and has a maximum when the satellite approaches the 



perigalacticon. The motion of the LMC through the hot halo 
of the MW during the last 1 Gyr is modeled using "test wind 
tunnel" simulations with increasing ram-pressure values. 

The paper is structured as follows. Section 2 describes 
the models and the star formation criteria adopted. Section 3 
illustrates the results of simulations without star formation, 
focussed on the investigation of pure effects of compression 
on the LMC interstellar medium while Section 4 describes 
the runs where star formation is activated. Several simula- 
tions have been performed, assuming different star forma- 
tion models, disk inclinations and hot halo densities. 



2 SIMULATIONS 
2.1 Galaxy model 

The initial conditions of the simu lations are constructed 
using the technique described by iHernauistI (|l993l ). Our 
disk galaxy model is a multi-component system with a 
stellar and gaseous disk emb edded in a spherical NFW 
(jNavarro. Frenk Whitelll99/t ) dark matter halo. The den- 
sity profile of the dark matter halo is adiabatically con- 
tract ed in response to baryonic infall (|Blumenthal et aP 
1986). The stellar disk follows an exponential surface density 
profile of the form: 

m) = ,z^e^p(-R/R^), (1) 

ZTvKd 

where Md and Rd are the disk mass and radial scale length 
(in cylindrical coordinates), respectively, while the thin ver- 
tical structure has a scale length Zd ^ ^Rd- 

p^{R^z) = ^^sech'{z/zd). (2) 

^'Zd 

The gaseous disk is characterized by an exponential profile 
with the same radial and vertical scale length as the stellar 
component and by a constant density layer which extends 
up to SRd. 

The structural parameters of the disk and the halo 
are chosen so that the resulting rotation curve resembles 
that of a typical bulgeless late- type ( Sc/Sd) disk galaxy 
(|Courteaulll997l : IPersic Saluccil [l997h . They are similar 
to those adopted in M05 for the initial LMC model and re- 
produce quite well the pe ak of the rotation curve inferred by 
van der Marel et al.1 (|2002 ) (Fig. 1). As seen in M05, the in- 
teraction with the MW does not affect significantly the stel- 
lar and dark matter mass in the inner 8 — 9 kpc of the LMC 
and consequently the global rotation curve within this radial 
range. The choice of an extended gaseous component for the 
initial LMC model is motivated by the fact that spiral galax- 
ies in the local universe are commonly observed to be em- 
bedded in extended disks of neutral hydrogen significantly 
larger than their stellar component (Hunter & Galla gher! 
Il985l : iBroeils k van WoerdenI Il994l ). As seen in M05, the 
combined effect of tidal interactions and ram-pressure strip- 
ping can remove a significant fraction of gas from a LMC 
disk orbiting within the hot halo of the MW, with a ram- 
pressure stripping radius which is a factor of three smaller 
than the initial radius of the gaseous disk. Also in the 
case of a LM C with orbital velocities significantly higher 
(Mast roDietrd l2008) hydrodynamic and gravitational forces 
together are effective in resizing and reshaping the extended 
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gaseous disk of the satellite beyond 8 kpc. In the present 
work we neglect the presence of gravitational forces focus- 
ing on the effects of pure ram- pressure. Therefore we do not 
expect to see a significant decrease in the radius of the gas 
distribution. However, in order to take in account the loss 
of cold gas from the disk of the satellite and the star forma- 
tion events, we assumed an initial amount of gas in the disk 
which is about 3 times larger than the HI mass in the LMC 
(2.9 X 10^ Mq according to Putman et al. 2003). 

The mass within the virial radius is set equal to 2.18 x 
10^° M0 and the fraction of mass in the disk is ^ 10%, 
equally distributed between the gaseous and stellar compo- 
nent. The contribution of the different components to the 
global rotation curve, assuming a disk scale length Rd = 1.7 
kpc and a dark halo concentration c = 9.5 (where c is de- 
fined as c = Vvir/vs, with rvir and rs virial and scale radius 
of the NFW halo, respectively) is plotted in Fig.[T] The halo 
spin parameter, which sets the disk scale length in our mod- 
eling, is A = 0.074, where A relates the angular momentum 
J and the total energy ^ of a system with virial mass M^^r 
through the relation A = J\E\^^'^G-^ M'^^'^ 

The initial stellar disk of the satellite galaxy has, within 
its scale radius Rd, a central mass surface density of ^ 
35M0pc~^ (Fig- ^1 that corresponds to a B-band sur- 
face brightness of ^ 24 mag arcsec"^, assuming a mass to 
light ratio c^i 2. The central gas surface density is only 
^ 16M0pc~^ since a significant fraction of gas is dis- 
tributed in the external disk. Assuming 72% HI abundance 
this value corresponds to an hydrogen column density of 
~ 1.5 X 10^^ within Rd, compar able with the values observed 
bv lStavelev-Smith et al.] (|2003h with the LMC Parkes multi- 
beam HI survey. 

2.2 Stability criterium 

In order to obtain a strongly stable disk against bar for- 
mation even in the presence of significant gas stripping and 
consequent perturbation of the satellite potential, the thick- 
ness of t he stellar component is set such that the Toomre's 
(|Toomra 1964) stability criterion is largely satisfied. In par- 
ticular the Toomre's parameter for the stellar disk: 



Qs{R) 



ar{R)K{R) 

3.36GS,(i?)' 



(3) 



where ar{R) is the radial velocity dispersion, k{R) is the 
local epicyclic frequency and Ss(i?) the unperturbed stellar 
surface density, has a minimum at the disk scale length with 
Qs{Rd) ~ 4. For a gaseous disk the stability of the disk is 
expressed in terms of the gas sound speed Vs and surface 
density S^(r) through the relation: 



Qg{R) 



VsK{R) 



(4) 



7rGi:g{R)' 

The gaseous disk has initially a constant temperature of 

10000 K, which imp lies Qa(R) > 3 a nd Q(Rd) = 3.2^ 

According to Ijog SolomonI |l984) and iRafikovl 

(|200lh . the stability of a multicomponent disk is not guar- 
anteed by the individual stability of its single constituents, 
due to the mutual gravitational interaction between gas 
and stars. Stars are characterized by velocity dispersions 
3-4 times larger than the typical sound velocities in the 
cold gaseous disk and even relatively small variations of the 




Figure 1. Galaxy model rotation curves. 

gaseous component paramet ers can significantly a ffect the 
stability of the whole disk (I Jog SolomonI 1 19841 ). There- 
fore we choose a large value of Q to contrast the effects of 
ram-pressure. In the case of a two components - gaseous 
and stellar - disk, th e stability condition is expressed by 
(Ijog Solomonlll984l ) 



Qt 



3.36 1 + ^2 1 + ^2^2 



> 1, 



(5) 



where q = kas/i^ and x — Cg/ag. Fig. [3] illustrates the de- 
pendence of Qtot on the dimensionless wavenumber of the 
perturbation q within three different regions of the disk: at 
the disk scale length Rd^ R — kpc and in the external 
regi on (j^ = 8 kpc), where the gas component predominates. 
The Jog & Solomon (1984') criterium is always satisfied and 
the disk is stable against axisymmetric perturbations, inde- 
pendently of their wavelength. In order to check stability the 
disk was initially evolved in isolation for 1 Gyr. 



2.3 Star formation recipes 

All the simulations we now discuss were carried out us- 
ing GASOLIN E, a parallel tre e-code with multi-stepping 
{ Wadslev, Sta del Quinnl2004l ) which is an extension of the 
pure N-body gravity code PKDGRAV developed by IStadel 
(I2OOII). The code uses a spline ker nel with compact sup- 
port (jMonaghan Lattanzij Il985l l where the interaction 
distance for a particle i is set equal to two times the smooth- 
ing length hi^ defined as the k-lh neighbour distance from 
the particle. In this paper A; = 32. The internal energy of 
the gas is inte grated using the antisymmetric formulation of 
lEvrardl (|l988l l that conserves entropy closely. Dissipation in 
shock s is modeled using the quadratic term of the standard 
Mon aghanI (j 19921) artificial viscosity. The Balsara (Balsaraj 
I1995I ) correction term is used to suppress the viscosity in 
non-shocking, shearing environments. The code includes ra- 
diative cooling for a primordial mixture of hydrogen and he- 
lium in collisional equilibrium. At temperatures below 10^ 
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Figure 2. Gaseous and stellar disk surface density profiles (cylin- 
drical coordinates). 
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Figure 3. Total stability parameter Qtot of the disk (see Equa- 
tion [5]) versus the dimensionless wavenumber q = kas /^c of the 
perturbation. Qtot is calculated for three different radii of the 
disk. 



star formation to collapsed, virialized regions. The physical 
density threshold describes the steep drop in star formation 
rate observed in disk galaxies when the gas surface density 
is much lower than a critical value T>c (Kennicutt 1989). The 
density threshold pc = 0.1 c ni~^ is compatible with obser- 
vational results. According to iKat j (|l992l ). the star forming 
region has to be part of a converging flow that implies a local 
negative divergence of the SPH velocity field. However, the 
converging flow criterium was introduced to describe star 
formation in cosmological simulations, where the geometry 
of the collapsing regions is approximatively spherical. In the 
case of star formation regions like 30 Doradus, localized at 
the periphery of the LMC disk where gas particles relatively 
close in distance can have significantly different kinematics, 
this criterium leads to underestimate of the star formation 
rate. Therefore the converging flow is not required in most of 
the simulations. A single run including the converging flow 
requirement has been performed for comparison. 

The star formation rate is assumed to be propor- 
tional to pgi^s (|Silkl [l987l ). where pgas represents the vol- 
ume density of the cold gas, and is given by the expression 
(Katz, Weinbe rg Hernquist..l996. ) 



dp^ 
dt 



C-kPgas 
tform 



(6) 



where the star formation timescale tform is the maximum 
between the local gas dynamical collapse time tdyn — 
{4:7TGpgas)~^^'^ and the local cooling time. If the gas is al- 
ready cool enough to form stars i.e T < Tmax, then tdyn is 
used. We assumed Tmax = 30000K. The constant star for- 
mation rate parameter is chosen such that we reproduce 
the global LMC star formation rate (Sandage 1986). Once 
a gas particle satisfies the above criteria, it spawns stars ac- 
cording to a probability distribution function. In particular, 
the probability p that a gas particle forms stars in a time 
At is modeled as 



(1 



(7) 



A random number is then drawn to determine whether the 
gas particle forms stars during At. For all the simulations 
in this paper At = 4 Myr. The newly created collisionless 
particle has the same position, velocity and softening length 
as the original gas particle while its mass is a fixed fraction e 
of the parent gas particle, whose mass is reduced accordingly. 
Following Katz (1992) we assumed for our favorite models a 
star formation efficiency e — 1/3. Up to six particles are then 
created for each gas particle in the disk. After its mass has 
decreased below 10% of its initial value the gas particle is 
removed and its mass is re-allocated among the neighbouring 
particles. 



K the gas is entirely neutral and due to the lack of molecu- 
lar cooling and metals, the efficiency of the cooling function 
drops rapidly to zero. We used a star formation recipe that 
includes density and temperature criteria, while converging 
flow criterium is not required in most of the simulations. 
Gas particles are eligible to form stars only if the density 
of the star formation region has a minimum phy sical den- 
sity corresponding to 0.1 hydrogen a toms per cm^ ([Summersl 
ll993l : lKatz. Weinberg Hernauist ll996) and an overdensity 
Pgas/ pgas > 55.7 (Katz ct al. 1996), which basically restricts 



2.4 Test wind tunnels 

In order to study the influence of pure ram pressure on a 
galaxy model orbiting in a Milky Way halo, we performed 
"wind tunnel" simulations where the ram-pressure value 
varies with time. 

We represent the hot g as as a flux of particles moving 
along the major axis of an oblong of base equal to the di- 
ameter of the dark matter halo of the satellite and height 
h — vt^ where v is the velocity of the LMC at the peri- 
galacticon and t is the time scale of the simulation. The hot 
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particles have an initial random distribution and a temper- 
ature T = 10^ K. The box has periodic boundary conditions 
in order to restore the flow of hot gas that leaves the oblong. 
The galaxy model is at rest at the center of the o blong. 

Kinematical data (|van der Marel et al.l I2OO2I I indicate 
that the LMC, presently located at ~ 50 kpc from the Galac- 
tic center, is just past a perigalactic passage and has an or- 
bital velocity of about 300 km s~^ . Recent proper motion 
measurements bv iKallivavalil et al. I (|2Q06 ) and Piatek et al. 
(|2007l l suggest that the velocity of the satellite is substan- 
tially higher (almost 100 km s~^) than previously estimated 
and consistent with the hy pothesis of a first passage about 
the MW (B esla et al.ll2007l ). In both scenarios the Cloud is 
affected by the largest ram-pressure values during the last 
million years of its orbital evolution. I ndeed, while in the 
models proposed by "Besla et ajj (|2007l ) the LMC does not 
enter the halo of the MW earlier than 1 Gyr ago (slightly 
different orbits are found in Mastropietro 2008), in M05 we 
have shown that the change in the orbital parameters due 
to dynamical friction strongly affects the ram-pressure strip- 
ping rate. Even in the case of a "low velocity" model the 
largest ram-pressure on the satellite is expected during the 
last orbital semi-period (about 1 Gyr) due to the increasing 
velocity and external gas density. 

We followed ram-pressure acting on the LMC's IGM 
during the past 1 Gyr. The density of the hot external gas 
increases with time, in such a way that in our low veloc- 
ity model the external pressure experienced by the cold 
disk varies from Pmin — pminV^in = 5 x 10~-^^dyncm~^ 
to Pmax — 1.52 X 10~^^dyncm~^ at the time correspond- 
ing to the pericentric passage. This is equivalent of as- 
suming Vmin — 170kms~^ and Vmax 300 km s""*^, and 
a number density of the external gas that increases from 
~ 10~^cm~^ to 10~^cm~^ at ~ 50 kpc from the Galactic 
center. These density values are comparable with those pro- 
vided by M05, who modeled the MW hot halo assuming 
a spherical distribution of gas that traces the dark matter 
profile, with a mean number density of 2 x 10~^ within 150 
kpc. We also consider the eventuality of a less dense Galactic 
halo and performed runs where the gas d ensity is a facto r 
10 lo wer. Models with higher velocities (jKallivavalil et al] 
20061) and o r bital parameters similar to those suggested by 
Besla et al.1 (|2007l ) are also explored. In details, Pmin is 
the same as in the low velocity models since the higher 
orbital velocity at the beginning of the simulation (about 
250 km s~^) is compensated by a lower external density 
(^ 5 X 10~^cm~^, according to Mastropietro 2008. Indeed 1 
Gyr ago the satellite has just passed through the virial ra- 
dius of the MW) . The maximum pressure felt by the disk is 
Pmax = 2.67 X 10~^^dyncm~^, that corresponds to a Cloud 
moving with Vmax 400kms~^ through an external hot 
medium of density 10~^cm~^. 

Each galaxy model is simulated using 750000 particles, 
of which 6 X 10^ are in the dark halo and 1.5 x 10^ in the 
disk (10^ collisional and 5 x 10^ collisionless) . The hot gas in 
the "wind tunnel" has 2 x 10^ particles, in such a way that 
the mass ratio mh/rridisk between hot particles and particles 
in the disk is close to the unity even when the halo density 
is the largest. This choice permits to avoid the presence of 
scattering and numerical holes which artificially change the 
shape of the front edge and influence the morphology of the 
disk (M05). The gravitational spline softening is set equal 



Figure 4. Density (upper panel) and temperature (bottom) map 
of a 50x50 kpc^ region around the satellite, which is moving 
towards the left of the page nearly face-on (inclination angle 
i = 10°, cooUG). Logarithmic color scaling is indicated by the 
key at the bottom of the figure. The surface density map ranges 
between 2 x 10^ (black) and ~ 2 x 10^ M© kpc~2 (yellow), 
while temperatures vary from 10"^ (black) to 10^ K (yellow). 



to 0.5 kpc for the dark halo and the hot gas in the oblong, 
while it is 0.2 kpc for stars and gas in the disk. 



3 COOLING SIMULATIONS 

In order to study the effect of pure compression on the den- 
sity distribution of cold gas in the LMC disk, we have run 
a first set of simulations where the gas cools radiatively but 

star formation was not activated. 

van der Marel et all (I2OO2I I and 
angle between 



According to 
IKallivavalil et al.1 (|2( 



the present 

the LMC's disk and its proper motion vector is roughly 
30°. Even neglecting the effects of precession and nutation 
on the disk plane of the satellite (^van der Marel eralll2002l ) 
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Figure 5. Gas temperature (top panel), density (middle) and 
one dimensional velocity (bottom) profiles across the shock dis- 
continuity. 

this angle is expected to vary significantly during an orbital 
period, especially in the proximity of a pericentric passage 
due to rapid changes in the velocity vector. DiflFerent 
relative orientations of the disk with respect to the orbital 
motion are therefore investigated. The inclination angle i is 
defined as the angle between the angular momentum vector 
of the disk and the flux of hot particles in the wind tube, so 
that a galaxy moving edge-on through the external medium 
is characterized by i = 90°, while the observed LMC disk 
would have i ^ 60°. We explored cases with inclination 
angle i of 90, 45 and 10° (runs cool90, cool45 and cooUO, 
respectively) . 

With a hot halo temperature of 10^ K the relative ve- 
locity between the satellite and the external medium is su- 
personic (sound speed ~ 135 km s~^ and Mach number 
A4 — 2.2 and 3 at the pericenter of the low and high ve- 
locity orbit, respectively) and a bow shock forms in front of 
the disk (Fig. 2]). Since the cooling time of the post-shock gas 
is 16 Gyr, the shock can be considered adiabatic and hy- 
drodynamical quantities at the two sides of the shock front 
are in first approximation related by the Rankine-Hugoniot 
equations for a stationary normal shock. 

For A4 > 1 the jump conditions give pi/p2 —v^jvx < 1 
and Ti < T2 , where subscripts 1 and 2 denote upstream and 
downstream quantities. The ram-pressure P2 = P2V2 actu- 
ally felt by the galaxy behind the shock front is therefore 
smaller than that it would suflFer due to the upstream flux 
of hot particles, but conservation of momentum flux across 
the shock discontinuity implies that the reduction in dynam- 
ical pressure has to be balanced by an increase in thermal 
pressure (see also Rasmussen et al. 2006). 

Fig. [5] illustrates the behavior of hydrodynamical quan- 
tities across the shock discontinuity for a snapshot corre- 
sponding to the perigalacticon of a low velocity orbit. The 
disk inclination is i = 10°. The horizontal axis is centered 
on the LMC stellar disk and oriented perpendicularly to 



Figure 6. Dynamical (Pdyn), thermal (Ptherm) and total (Ptot) 
pressure across the shock. 

the bow shock nose, with the shock located at x 10 kpc 
and the satellite moving towards increasing values of x. The 
Mach number derived by the temperature jump is M = 2.1, 
in good agreement with the theoretical value for a normal 
shock. Only hot halo particles are considered in computa- 
tions but, due to the SPH nature of the simulations, close to 
the border of the disk we observe a further density increase 
and a sharp drop in temperature. The x- velocity proflle is 
plotted in the system of reference where the pre-shock gas 
is at rest. 

Pressure profiles across the shock are plotted in Fig. [6l 
where Pdyn is defined as pv^. The total pressure remains 
roughly constant until the edge of the disk. The steep in- 
crease at X = 2.5 kpc is due the rapid growth in density at 
the border of the disk, which is not immediately followed by 
a decrease in temperature. 

For i > 0° the shock wave is inclined with respect to 
the initial flow velocity and the Rankine-Hugoniot condi- 
tions apply to the normal components of the velocity across 
the shock discontinuity, while the component parallel to the 
shock front remain unchanged. The flow is therefore de- 
flected toward an oblique shock wave and the jump at the 
shock discontinuity is smaller. 

An edge-on {i = 90°) disk behaves like a wedge moving 
supersonically with the vertex facing upstream. If the wedge 
angle is smaller than or equal to the maximum flow deflec- 
tion angle, the oblique shock becomes attached to the vertex 
of the wedge and the flow is deflected so that the streamlines 
are parallel to the surfaces of the wedge. The shock stand- 
ofl[ distance thus depends on the external density proflle of 
the collisional edge-on disk and on the Mach number of the 
incident flow. In both low and high velocity edge-on models 
the shock results almost attached to the disk. 

Figs. EHSl illustrate - for diflFerent values of the inclina- 
tion angle i - changes in the disk gas density distribution as 
the satellite passes through increasing values of the external 
pressure, moving towards the perigalacticon. Each couple of 
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Figure 7. Edge-on model (run cool90). Evolution of the disk gas density distribution. Left panels: HI column density maps (the galaxy 
is moving towards the left of the page). The color-scale is logarithmic with limits 5 x 10-*^^ and 1.2 x lO^-*^ cm~^. Right panels: mean gas 
density and radius of the external disk as a function of the azimuthal angle cf). Each pair of plots (from the top to the bottom) represents 
the state of the disk at increasing times along the orbit: Time= 0.4, 0.6, 0.8, IGyr. 
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Figure 8. The same as in Fig. [71 for the model with inchnation angle i = 45° (run coo 
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Figure 9. The same as in Fig. [Tlfor the model with inchnation angle i = 10° (run cooUO). 
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panels illustrates the state of the disk at increasing (from 
the top to the bottom) times along the orbit. Panels on the 
left represent HI column density maps. The density contrast 
is chosen in order to highlight the density gradient in the ex- 
ternal disk, since the gas distribution in the central regions 
of the LMC is dominated by the presence of the bar and 
a direct comparison with pure ram-pressure simulations is 
not possible. The color scale is logarithmic, with white cor- 
responding to a density larger than 1.2 x 10^^ cm~^ and blue 
to values lower than 5 x 10^^ cm~^. Hot gas particles flow 
on to the disk from the left to the right of each plot, with 
increasing ram-pressure values from the first to the fourth 
image. The disk is seen face on and rotating clockwise. In 
the case of the edge-on run (Fig. [Tj) cold gas particles lying 
in the left-bottom quadrant of the disk feel the largest ram- 
pressure, due to the fact that their relative velocity with 
respect to the external medium is maximum. The rotational 
velocity of the external disk is ^ 55 km s~^, which implies 
a relative velocity at the pericenter of ~ 350 km s~^ (450 
km s~^ in the case of a high velocity orbit). Panels on the 
right represent the change in mean density and radius of 
the gaseous external disk as a function of the azimuthal an- 
gle (/). Referring to the geometry of the HI density images 
on the left, = 0° corresponds to the bottom of the disk 
and increases clockwise in such a way that the disk moves 
in the direction oi (f) — 90°. The gas density and the mean 
radius are both calculated within sections of a three dimen- 
sional annulus with internal and external radius equal to 7 
and 15 kpc, respectively. The initial azimuthal profiles (not 
represented in the plots) are flat since both these quanti- 
ties have only radial dependence. As soon as the satellite 
starts moving through the surrounding medium, the exter- 
nal gas density develops a peak centered on = 90°: disk 
particles localized in regions of maximum ram-pressure get 
compressed and move on inner orbits, while their circular 
velocity increases consequently. After about a quarter of the 
orbital period the gas has reached its minimum radius and 
maximum local density. The gaseous disk becomes strongly 
asymmetric: compression at the front edge produces a den- 
sity increase along the left border of the disk, evident in 
the HI maps even at early times. The high density region 
forms a thin 1.5 kpc) but continuous and well defined 
arc which has not an equivalent in the stellar distribution. 
At the perigalacticon this feature extends for almost 160° 
with a density more than one order of magnitude higher 
than gas located at smaller radii. Its average thickness (^ 
1.5 kpc) and velocity dispersion along the line of sight are 
larger than the average values in the rest of the disk. 

In the case of a satellite moving through the hot medium 
with an inclination angle different from 90° the external 
pressure directed perpendicularly to the p lane of the disk 
mcreases as cost (Roediger k Briiggenl lioOGl ) while compres- 
sion at the leading edge is much less pronounced. Figs.[8]and 
[9] refer to runs with inclination angles i — 45° and i — 10°. 
Disks are shown face-on. The increment in density along the 
leading edge is smaller (cool45) than in the edge-on model 
and almost absent for i — 10° (coollO), while compression 
perpendicular to the plane of the disk produces local gravita- 
tional instabilities in the external gaseous disk (also Mayer, 
Mastropietro & Tran in preparation). This effect is more 
evident in the nearly face-on run coollO where high density 
filaments delimitate regions where the local density is almost 
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Figure 10. HI column density profile of the final LMC disk in 
units of lO^'-'cm"^. The curves refer to different values of the in- 
clination angle i. The red long-dashed line represents the column 
density profile of the initial disk. 



one order of magnitude lower. Despite the absence of a peak 
in the azimuthal mean density profile, the integrated final 
density of coollO is comparable to the other runs. 

Fig. 1101 represents the azimuthally averaged HI column 
density profile of the final disk configuration for the three 
simulations. As a result of the increase in density along the 
edge of the disk, the mean column density shows a secondary 
peak at large radius. A li mb-brightened density profile has 
actually been observed by St avelev- Smith et al.l (|2003l l us- 
ing the Parkes multibeam HI survey of the LMC. The outer 
profiles of cool90 and cool45 are very similar, while in run 
coollO the increment in density with respect to the original 
profile of the disk (long-dashed red curve in the plot) is lo- 
cated at larger radii (^ 8.5 kpc). Indeed, due to compression 
of the leading edge the final gas distribution of cool90 and 
cool45 is asymmetric, the dense front edge being much closer 
to the center than the opposite border of the disk. There- 
fore the HI peak in the external disk is located at relatively 
small radii, while the azimuthally averaged gas distribution 
is more extended than in the case of the run coollO since 
ram pressure elongates the back side of the disk. Again we 
stress the fact that our gaseous disk is more extended than 
it would be in a fully self consistent simulation including 
both gra vitational and hydrodynamical forces. As seen in 
M05 and i Mastropietrol ((2008) the combined effect of ram- 
pressure and tidal stripping is quite efficient in stripping 
gas from the outer satellite's disk, creating the tip of the 
Magellanic Stream already at large distances from the MW. 
The high density feature in cool90 would then form along 
the border of the disk at smaller radii and would not be 
easily subjected to further stripping due to the relatively 
ram- pressure values. 
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Table 1. Star formation simulations. The second column indicates the use of the converging flow criterium in the star formation recipe, 
the third and the fourth columns represent the star formation efliciency e and the star formation rate parameter (Eq. (G]), respectively, 
while i is the inclination angle of the disk, i.e. the angle between the angular momentum of the disk and the flux of hot particles. A 
model with i = 90 is indeed a disk moving edge-on in the external medium. The sixth and seventh columns represent the minimum 
and maximum ram pressure suffered by the satellite during its approach to perigalacticon. The last column indicates the presence of a 
temperature floor for the disk gas. 
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4 SIMULATIONS WITH STAR FORMATION 

The compressive increase in HI density is naturally asso- 
ciated with excess star formation. SPH simulations cannot 
follow the formation of molecular clouds but in first approx- 
imation the molecula r g:as fraction can b e related to the 
density of atomic gas (|Vollmer et aljlioosl v 

The main parameters of star formation simulations are 
summarized in Tabled] 

As we already mentioned in Section 3, our standard star 
formation model (SF) does not include the converging flow 
criter ium and is characterized by an efficiency e = 0.33 ( K at z 
I1992I V The star formation rate parameter is initially set 
equal to 0.02. This model was adopted to run wind tube sim- 
ulations with inclination angles i = 90, 45, 10° (SF90, SF45, 
SFIO) (Table [T|). We also investigated different star forma- 
tion recipes requiring converging flows and assuming differ- 
ent values of and e. We explored star formation rate pa- 
rameter values in the range from 0.01 to 0.05, that produce 
a star formation rate integrated over th e entire d isk compa- 
rable with the O.IM0 yr"^ provided bv lSandagj (jl986). An 
efficiency e — 1 implies that whenever a gas particle satisfies 
the star formation requirements, it is im mediately turned 
into a single star particle of the same mass (|Kaufmann et all 
I2OO6I ). 

In the last six simulations listed in Table [1] we used our 
standard star formation model SF to investigate the effects 
of different orbital parameters and gas halo densities. 

Runs SFv400 are characterized by a maximum ram 
pressure value corresponding to a perigalactic velocity of 400 
km s~^ (|Kallivavalil et al]|2006h . Such high velocity disks - 
when moving face-on through the external hot medium - 
are strongly affected by local instabilities and star forma- 
tion is consequently enhanced. The introduction of an arti- 
ficial lower limit for the satellite gas temperature (in runs 
SF10v400t 12000 and SF10v400t 15000), higher than the cut- 
off in the cooling function, has the effect of reducing grav- 



itational instabilities and fragmentation in the disk. This 
temperature threshold can be justified in order to crudely 
model th e effect of the UV background and stellar feedback 
(|Barnesl[20Q2. ). 

Finally, with simulations SF901d and SFlOld in Table □ 
we also consider the possibility of a Galactic hot halo ten 
times less dense than our standard model. 

Fig. [11] illustrates the state of the newly formed stellar 
disk at increasing times along the satellite orbit. Each couple 
of rows corresponds to one of the first three runs of Table 
[T]and is associated with a different inclination angle i. The 
first row of the pair represents the face-on projection of the 
disk with the galaxy moving towards the left of the page 
and the same geometry as in Fig. [7] Time increases from 
the left to the right. Each small cross indicates a new star 
formation event at the time of the snapshot (within a time 
interval of 40 Myr) while the circle delimitates the external 
disk (r > 7 kpc). Stars form in the central regions as soon 
as the star formation algorithm is activated, but here the 
star formation activity of the inner disk is not represented. 
The second row represents the total mass of the newly 
formed stars in the external disk as a function of the disk 
azimuthal angle 0. 

In the case of a galaxy moving edge-on through the ex- 
ternal medium (SF90), stars form at the leading edge of the 
disk when the ram-pressure becomes larger than ~ 5 x 10~^^ 
dyn cm~^, at Time ~ 0.3 Gyr. The location of the star for- 
mation events initially corresponds to the HI column density 
peak observed in cool90 around 100° (Fig. [7]). Later on it 
expands along the entire front edge, creating a thin stellar 
arc well distinct from the star formation events that char- 
acterize the central disk. As soon as the satellite encounters 
ram pressure values comparable to those experienced by the 
LMC at the perigalacticon (Time >0.6 Gyr) some episodes 
of star formation occur even on the back side of the disk 
(last plot on the top right of Fig. [TT]) , although they are 
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Figure 11. Location of the instantaneous star formation events in the external disk. Each couple of rows corresponds to a different 
inclination angle of the disk. Time increases from the left to the right: Time= 0.4, 0.6, 0.8, 1 Gyr. The first row of the couple represents 
the disk face on projection, the second one the star formation rate versus the azimuthal angle 0. 
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not relevant in terms of new stellar mass formed. Indeed 
shows a drastic drop at = 200°. 

Runs with inclination i < 90° are characterized by 
significant star formation only for values of the external 
pressure larger than 10~^^ dyn cm~^ . In the case of the 
nearly face-on run SFIO, at t > 0.7 Gyr star formation 
occurs in the entire external disk. Ram-pressure affects 
the plane of the disk almost perpendicularly and stars 
form along the delocalized and filamentary high density 
structures visible in Fig. [9] (M ayer et al. in prep . ). Con trary 
to what has been found by ..Kronberger et al.l (|2008h who 
focused on higher ram-pressure values (n ~ 10~^ cm~^ and 
V = 1000 km s~^) typical of the outskirts of galaxy clusters, 
the newly formed stars are all located in the plane of the 
satellite's disk (with the exception of the high velocity 
face-on run SF10v400 where about 10% of the stars forms 
behind the disk). The star formation events appear to 
be distributed nearly homogeneously along the azimuthal 
profile of the external disk, although a small peak in is 
observable near i — 90°. In fact, the orientation of the disk 
with respect to its orbital motion is not exactly face-on. 
The case of the intermediate run SF45 is more complex. In 
a first phase, for low ram pressure values, star formation 
is produced by compression at the leading edge and a thin 
star formation front - although not so well defined as in the 
case of a pure edge-on model - appears on the east side of 
the disk. As soon as the external pressure reaches a critical 
level compression directed perpendicularly to the disk 
becomes the dominant mechanism driving star formation. 

Converting SFRs to Ha luminosities according to 
iKennicuttI (|l998l ): 



L{Ha) 



SFR(M0yr-^) 
1.26 X 1041 ' 



(8) 



where SFR is the star formation rate averaged over the 
last 40 Myr (nearly two times the stellar age of 30 Doradus), 
we obtain the H-a maps illustrated in Fig. 1121 

High emission regions are mainly concentrated in the 
external disk (with the exception of run SFld90 where the 
ram-pressure exerted by the low density halo is not enough 
to induce star formation at the edge of the disk). The con- 
tinuous stellar arc forming along the leading side of the 
disk in edge-on runs breaks up into several distinct and 
very luminous H-a regions that more closely resemble the 
star-forming complexes observed on the eastern border of 
the LMC. The inclusion of stellar and supernovae heating - 
which has been neglected in the present simulations - could 
prevent further star- format ion around highly emitting re- 
gions and consequently produce more compact and isolated 
star-forming complexes. Nevertheless, modeling single star- 
formation complexes - whose linear extension is smaller than 
our softening length - is beyond the scope of this paper. 

The present inclination of the LMC's disk with respect 
to the orbital motion is about 60° (according to the conven- 
tion adopted in this paper). Since the satellite is currently 
near a perigalactic passage we expect the H-a map at the 
leading border of the disk to be something in between pure 
edge-on runs and the run with inclination of 45°. On the 
other side, it is very likely that the disk inclination during 
the ph ase of approach to t he pericenter was different. In- 
deed in iMastropietrol (|2008l l we have simulated the LMC's 
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Figure 13. Star formation rate of the external disk cf) for the 
first ten runs listed in Tabled Time increases approaching to the 
present position. 



orbit according to the n ew proper motion measurements of 
iKallivavalil et aP (|2006h and found that the cloud enters the 
MW halo face-on and moves almost face-on during most of 
the last 1 Gyr. It turns nearly edge-on only at the peri- 
galacticon. This would have a remarkable effect on the star- 
formation history of the external disk during the last 1 Gyr 
and some impact also on the H-a maps. Indeed, although 
the H-a emission would be mostly concentrated on the east- 
ern side of the disk due to the very recent edge-on motion, 
we expect to see some luminous clumps forming a patchy 
distribution on the entire disk, due to gravitational insta- 
bilities and subsequent star-formation induced by a nearly 
face-on compression of the gaseous disk before 30 Myr ago. 

The high velocity edge-on run SF90v400 presents a 
more elongated and thinner stellar arc along the leading bor- 
der, with a geometry similar to that obtained increasing the 
star formation rate parameter to 0.05 (SF90c0.05). The H-a 
map of SF90v400 (third panel of the third raw) shows two 
distinct luminosity peaks. One is located at the south-east 
region of the disk, roughly corresponding to the position of 
SODoradus and the two compact emission regions N159 and 
N160. 

Fig. [13] illustrates how the star formation rate of the 
external disk changes with time in the different models. The 
three black curves refer to the standard star formation runs 
SF, characterized by low orbital velocities. The edge-on disk 
starts forming stars earlier, but for large ram pressure val- 
ues the star formation rate of SFIO and SF45 grows faster. 
At the perigalacticon passage SF45 has indeed a higher star 
formation rate than SF90. In the case of an isolated LMC 
model star formation is almost absent for r > 7 kpc. The 
remaining curves in Fig. [13] refer to different star forma- 
tion recipes (rows 4-9 of Table [T]). The location of the star 
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Figure 12. H-a maps of the LMC's disk. The different panels represent (from top to bottom and from left to right) runs listed in Table 
□ SF90, SF45, SFIO, SFconv, SFel, SFconvel, SFQOcO.Ol, SFlOcO.Ol, SF90c0.05, SFlOcO.05, SF90v400, SF10v400, SF10v400tl2000, 
SF10v400tl5000, SFld90, SFldlO. The color-scale is logarithmic with hmits lO^i and 10^4 erg s^^ kpc-^. 



formation events in the external regions of the disk does 
not change significantly choosing different parameters in the 
star formation algorithm. The star formation rate parame- 
ter determines the amount of new stars forming but does 
not affect the minimum threshold in ram-pressure neither 
the evolution of the star formation rate. In particular, in 
the case of edge-on runs after an initial steep increment the 
curve seems to converge to a constant value for increasing 
external pressures. The consequences of an increased star 
formation efficiency (el) are almost negligible (but not in 
Ha maps where only the very recent star formation rate is 
taken in account: compare the first panels of the first and 
second row) while including the convergency requirement 



(SFconv) has nearly the same effect than reducing the star 
formation rate parameter of a factor two. 

For convenience, the star formation rate of the last six 
runs of Table [T] is plotted separately fFig. I14|) . The high ve- 
locity edge-on model SF90v400 is characterized by a steeper 
increment in star formation at earlier times but later on 
the curve flattens and the star formation rate at the peri- 
galacticon is similar to that of the low velocity case SF90. 
Differences in H-a maps are produced by a difference of 
5 X 10^ M© Myr"^ about Time = 1 Gyr. On the contrary, 
the star formation generated by compression perpendicular 
to the disk increases with increasing ram-pressure values (it 
shows a decrement only towards the end of the simulation) 
and in the case of the high velocity run SF10v400 reaches a 
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Figure 15. Star formation rate of the entire disk for selected runs 
in Table [T] SFiso refers to the LMC model run in isolation. 



peak ^1.5 times higher than in SFIO. The star formation 
rate of SF10v400 is strongly affected by the introduction 
of an artificial temperature threshold. With a temperature 
floor of 15000 K we nearly suppress star formation in the 
external disk. However, a threshold of 12000 K is already 
very high for a LMC model and more typical of luminous 
disk galaxies like the MW. A temperature floor lower than 
10000 K would not make sense since below this tempera- 
ture the cooling function adopted in th e present paper drop s 
rapidly to zero. If, as pointed out bv iMastropietrol (|2008h . 
the satellite is moving almost face-on until it gets very close 
to perigalacticon, we would not expect to see star formation 
before 0.6 Gyr independently of the orbital velocity. Differ- 
ences between high and low velocity runs should be marginal 
also near the pericenter since star formation in edge-on runs 
- assuming our standard prescriptions for star formation - 
seems to saturate around ~ 4 x 10^ Mq Myr~^. Finally, a 
hot halo ten times less dense than the one assumed in our fa- 
vorite model would reduce drastically the star formation in 
a face-on LMC and suppress completely the star formation 
on the leading edge of the disk. 

The star formation rate of the entire disk is illustrated 
in Fig. 1151 This plot is only indicative since we neglect the 
presence of the bar and its influence on the star formation 
history of the satellite. The star formation rate of a LMC 
model evolved in isolation is plotted for comparison. The 
largest contribution to SFR is given by star formation events 
in the central region of the disk. Curves peak between 15 and 
35 M© Myr~^, comparable with observations of the Mag- 
ellanic Clouds Photometric Survey (Harris et al. in prep). 
The total star formation rate is clearly not affected by ram- 
pressure before 0.4 Gyr. Indeed, after an initial sharp incre- 
ment curves are rather flat despite orientation and intensity 
of the external pressure. An initial burst in star formation is 
obtained only by increasing the star formation efficiency pa- 



rameter to 0.05. For Time> 0.5 Gyr the largest deviations 
from an average star formation rate of ^ 15 Mq Myr~^ 
are produced in the edge-on low density run SF901d, where 
dM/dt drops to 10 M© Myr"^ at Time = 1 Gyr, and in the 
high velocity face-on model SF10v400 whose star formation 
rate increases up to ^ 24 Mq Myr~^ toward the end of the 
simulation. Differences among the other models are of the 
order of few Mq Myr~^ and vary with time so that it would 
be quite difficult to use the star formation history of the en- 
tire disk to test the Cloud's orbital parameters and the hot 
halo density. 

Fig. [T7] represents the mean stellar age of the external 
disk versus the azimuthal angle for the same selected runs 
of Fig. 1151 The maximum increment in age in a clockwise 
direction is associated with edge-on runs, where stars form- 
ing at the leading edge move, in time, away to the side, due 
to the clockwise rotation of the disk. The youngest stars are 
located at 30° < < 100°. Clearly the gradient in age is 
much weaker in models with i < 90°. SF10v400 - with a 
mean stellar age of ^ 250 Myr - forms stars earlier with 
respect to the other nearly face-on runs, while the difference 
between SF90v400 and the corresponding low velocity run 
SF90 is about 5 Myr at the leading edge and not significant 
in the rest of the external disk. 

For the same runs we also plotted the final radial gas 
density profile (Fig. llSp . In most of the cases a secondary 
peak is still present. 



5 CONCLUSIONS 

We have performed high resolution "wind tunnel" simula- 
tions to study the effects of ram-pressure by a tenuous Galac- 
tic hot halo on the HI morphology of the LMC's disk, its re- 
cent star formation history, and location of the youngest star 
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Figure 16. Age distribution of the newly forming stars in the LMC's disk. The different panels represent (from top to bottom and 
from left to right) runs listed in Table □ : SF90, SF45, SFIO, SFconv, SFel, SFconvel, SFQOcO.Ol, SFlOcO.Ol, SF90c0.05, SFlOcO.05, 
SF90v400, SF10v400, SF10v400tl2000, SF10v400tl5000, SFld90, SFldlO. Magenta indicates stars older than 725 Myr. 



forming regions. We did not focus on the mass loss produced 
by ram-pressure stripping since this would also be affected 
by tidal interactions. For the same reason our galaxies do 
not form any bar and we actually start with a galaxy model 
very stable against bar formation so that the pure effects 
of external pressure are more clearly visible. Our LMC is a 
multi-component system with a spherical NFW halo, an ex- 
ponential stellar disk and a gaseous disk that extends up to 8 
times the stellar disk scale length. In each simulation the ex- 
ternal flux of hot particles increases with time as the satellite 
approaches the perigalacticon in such a way that the pres- 
sure experienced by the disk is consistent with the LMC's 
orbital velocity and an average hot halo density of 2 x 10~^ 
cm~^ within 150 kpc from the Galactic center. Low velocity 



runs are char acterized by a "classic" p ericentric velocity of 
300 km s"^ (|MastroDietro et al.1 liooa l while high velocity 
runs have velociti es compatible with the new proper motion 
measurements of iKallivavalil et al ] 12006). We expect the 
angle between the LMC's disk and its proper motion to vary 
significantly during the last billion years of orbital interac- 
tion. We have defined the inclination angle i as the angle be- 
tween the angular momentum vector of the disk and the flux 
of hot particles in the win d tube, s o that the observed LMC's 
disk would have i ^ 60°. Mastro pietrol (j2008) performed self 
consistent Nbody/SPH simulations of the interacting system 
MW/LMC adopting orbital constraints from the last LMC's 
proper motion measurements and found that the Cloud en- 
ters the MW halo face-on (i ~ 0°), moving nearly face-on for 
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Figure 17. External disk mean stellar age in function of the 
azimuthal angle (j) for selected runs in Table [T] 
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Figure 18. Azimuthally averaged HI column density profile of 
the final LMC disk for a selection of the runs listed in Fig. ^ 



most of the last billion years. It turns edge-on only about 30 
Myr ago. This means that the LMC is moving nearly edge- 
on close to the perigalactic passage, corresponding to the 
maximum ram-pressure valu es, consistently wit h the actual 
disk inclination measured bv iKallivavalil et al ] (2006). 

We have performed several simulations varying the in- 
clination angle of the disk, the star formation recipe and the 
intensity of the external pressure. We have shown that: 

• The compression of the leading border of an edge-on 
LMC disk can account for the high density HI region ob- 
served at the south east. In our simulations this high den- 
sity feature is well defined (with a mean density one order of 



magnitude higher than the surrounding gas) and localized 
within 1.5 kpc from the border of the disk. Its average thick- 
ness and velocity dispersion along the line of sight are larger 
than the average values in the rest of the disk. In cool90 it 
extends for almost 160° and c ould also explain the origin o f 
the spiral arm E described bv lStavelev- Smith et al. 
which does not have an equivalent in the stellar disk. 

• Compression directed perpendicularly to the disk (in 
runs with i < 90°) produces local instabilities in the gas 
distribution and a clumpy structure characterized by voids 
and high density filaments similar to those observed by the 
Parkes multibeam HI survey (see Fig. 3 of Staveley- Smith 
et al. 2003) . If the satellite was moving nearly face-on in the 
past - and according to iMastropietrol ()2008l l this is likely 
to happened during most of the LMC/MW orbital history 
- ram-pressure could be responsible for the general mottled 
appearance of the HI disk. 

• As a result of the increase in density along the edge 
of the disk the mean HI column density shows a secondary 
peak at large radius, in agreement with observations. 

• The compression of the satellite's IGM is naturally as- 
sociated with induced star formation activity. We focussed 
on the external regions of the disk since the central parts 
of the real LMC would be dominated by the bar. Edge-on 
disks start forming stars earlier, but for large ram-pressure 
values the star formation rate of runs with i < 90° grows 
much faster. The high velocity edge-on model SF90v400 is 
characterized by a steeper increment in star formation at 
earlier times but later on the curve flattens and the star for- 
mation rate at the perigalacticon is similar to that of the 
low velocity case SF90. On the other hand, the star forma- 
tion generated by a compression perpendicular to the LMC's 
disk increases with increasing ram-pressure values and in 
the case of the high velocity run SF10v400 reaches a peak 

1.5 times higher than in SFIO. If the satellite is moving 
almost face-on until it gets very close to the perigalacticon 
we would expect not to see star formation before 0.6 Gyr in- 
dependently of the orbital velocity. Differences between high 
and low velocity runs should be marginal also near the peri- 
center since star formation in edge-on runs - assuming our 
standard prescriptions - seems to saturate around 4 x 10^ 
Mq Myr~^ 

• In edge-on models the star formation of the external 
disk is characterized by a thin stellar arc along the leading 
border, well distinct from the star formation events in the 
central disk. If SFR is converted in Ka luminosities, this arc 
breaks in several distinct and very luminous Ha regions that 
more closely resemble the star forming complexes observed 
on the eastern border of the disk. Although the Ka emission 
is mostly concentrated on the eastern side as a consequence 
of the very recent edge-on motion, we expect to see some 
luminous clumps forming a patchy distribution on the entire 
disk, due to gravitational instabilities and subsequent star 
formation induced by a nearly face-on compression of the 
disk before 30 Myr ago. As observed by de Boe r et aP (|l998h 
stellar complexes on the leading edge show a progression in 
age in the clockwise direction, but a face-on compression in 
the recent past of the LMC would circumscribe this trend 
to the youngest stellar regions, with age < 30-40 Myr. 
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